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ABSTRACT 

We present and discuss Spitzer and near-infrared II2 observations of a new bi-polar protostellar 
outflow in the Rosette Molecular Cloud. The outflow is seen in all four IRAC bands and partially 
as diffuse emission in the MIPS 24 fim band. An embedded MIPS 24 /im source bisects the outflow 
and appears to be the driving source. This source is coincident with a dark patch seen in absorption 
in the 8 /im IRAC image. Spitzer IRAC color analysis of the shocked emission was performed from 
which thermal and column density maps of the outflow were constructed. Narrow-band near-infrared 
(NIR) images of the flow reveal H2 emission features coincident with the high temperature regions of 
the outflow. This outflow has now been given the designation MHO 1321 due to the detection of NIR 
H2 features. We use these data and maps to probe the physical conditions and structure of the flow. 
Subject headings: ISM: jets and outflows — ISM: individual objects (MHO 1321) — methods: data 
analysis 



L INTRODUCTION 

Outflows and jets from young stellar objects (YSOs) 
accompany the early stages of star formation. Outflows 
can manifest themselves as jets and knots of shocked 
material visible at optical and near-infrared wavelengths 
and also molecular emission observable at longer wave- 
lengths. The outflowing material plays a role in removing 
the excess angular momentum from the YSOs allowing 
them to evolve into stars. Outflows are able to trace 
the history of mass loss and accretion of their driving 
sources. Studying the structure and properties of these 
flows may provide clues to understanding the connec- 
tion betwe en jets and the associate d wide angle molec- 
ular flows (jReipurth fc Ballvl [2OOII ). Additionally, this 
outflowing material interacts with its surroundings and 
may affect its environment, possibly regulating further 
star formation and cluster evolution. The energy and 
momentum inputted by outflows may disrupt the sur- 
rounding ambient gas, contribute to the turbulence i n 
the cloud, and affect chemical processes (|BiII^l20?)l . 
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lYbarra &: Ladal (|2009D developed a technique to study 
the thermal structure of shocked H2 gas using color anal- 
ysis of observations from the Spitzer InfraRed Array 
Camera (IRAC). Given the vast amount of Spitzer data 
available, this technique can be used to survey large re- 
gions and simultaneously find and analyze shocked emis- 
sion. The IRAC color analysis enables the construction 
of temperature maps of the shocked gas which may in 
turn be used to probe the interaction of outflow with 
its surroundings. These maps may also be used to com- 
pare the properties of outflow with those of simulations 
allowing a better understanding of the physics involved 
and estimating the energy and momentum inputted by 
outflows into their environment. 

The Rosette Molecular Cloud (RMC) is a star form- 
ing region located at a distance of 1.6 kpc. Near- 
infrared imaging studies have revealed nine embed- 
ded clusters across the c loud (jPhelps fc Ladal 119971 : 
iRoman-Zufiiga et al.l I2008D . Outflow activity in the 
cloud has been revealed through th e [S I I] narrowband 
imaging survey oflYb arra fc PhelpsI (I200I and the ^^CO 
survey of iDent et al.l ( 20091) . In an analvsis of the Spitzer 
IRAC images of the Rosette Molecular Cloud, we have 
discovered a structure with the morphology of a bipo- 
lar outflow that is visible in the images from all four 
IRAC bands. This structure can be seen in the images 
published bv lPoulton et al.l ()2008() although it is not dis- 
cussed in their paper. 

In this study we analyze the outflow using near infrared 
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(NIR) narrowband imaging of the flow to confirm the 
presence of shocked gas inferred from analysis of the the 
IRAC images. We im prove the IRAC color analysis of 
lYbarra fc Ladal (|2009n and use it to create temperature 
and column density maps of the outflow. Using both the 
NIR and IRAC data, we probe the physical conditions 
and structure of the outflow. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer IRAC and MIPS data reduction 

We used MIPS 24^m and IRAC 3.6-8.0 urn data 
from program 3391 (PI: Bonnel) available in the Spitzer 
archive. The IRAC frames were processed using the 
Spitzer Science Center (SSC) IRAC Pipeline vl4.0, 
and mosaics were created from the basic calibrated 
data (BCD) frames usi ng a custom IDL program (see 
iGutermuth et all (|2008[ ) for details). The MIPS frames 
were processed using the MIPS Data Analysis Tool 
{Gordon et al. 2005). 

2.2. Near-Infrared H2 observations and data reduction 

Near-Infrared, narrow-band observations of the out- 
flow were obtained with the Infrared Side Port-Imager 
(ISPI) on the Blanco 4 meter telescope at the Cerro 
Tololo Inter- American Observatory (CTIO). ISPI em- 
ploys a 2048 X 2048 HgCdTe Hawaii-2 array with a 10.25 
X 10.25 arcmin field of view and a plate scale of ^ 0.305" 
pixel"^. The outflow was imaged using the II2 1—0 S(2) 
2.034 ^m filter (A^ = 2.0336 ^m, AA/A = .007), H2 1-0 
S(l) 2.122 ^m filter (Ac = 2.1262 ^m, AA/A = .01), and 
Kcont filter centered at 2.1462 ^m. The telescope was 
dithered with a 20 point dither pattern with a integra- 
tion time of 60s at each dither position. The images were 
taken on the nights of 2008 December 17 and 2008 De- 
cember 19 with total integration time of 40 minutes in 
each filter. 

The raw images were flat fielded, corrected for bad pix- 
els, and linearized using the task osiris from the CTIO 
Infrared Reduction package. For each image, a sky frame 
was created by median combining the dithered images 
closest in time to the image. The IRAF tasks msctpeak 
and msc image, which are part of the IRAF Mosaic Data 
Reduction Package, mscred, were used to correct the im- 
ages for geometric distortions. The high order polyno- 
mial distortion terms were calculated with msctpeak us- 
ing the 2MASS Point Source Catalog as the reference 
catalog and the distortion correction was applied with 
msc image. The corrected images were aligned and then 
combined to form the final science images. 

3. RESULTS AND ANALYSIS 

Figure 1 shows the outflow in all four IRAC bands. 
The outflow appears as patchy regions of diffuse emission 
with an overall structure that is elongated and collimated 
in the E-W direction. Figure 2 shows the narrowband 
near-infrared emission images of the outflow. The NIR 
II2 knots coincide with the diffuse emission seen in the 
IRAC images. The NIR II2 images confirm the presence 
of shocked gas and the interpretation that this structure 
is an outflow. The eastern end appears to truncate at 
a bow shock. Slightly west of the bow shock, the H2 
images reveal a small scale chaotic structure followed by a 
more linear chain of knots. The western end of the outflow 



Table 1 

Positions and flux estimates for the NIR H2 knots of MHO 1321 



Knot R.A. (J2000) Dec (J2000) H2 1-0 S(2) H2 1-0 S(l) 

a 06:35:20.08 -1-03:56:37.9 4.1±2.0 7.6±2.2 

b 06:35:21.64 +03:56:28.8 2.2±1.4 8.9±2.2 

c 06:35:21.82 +03:56:32.3 3.9+1.9 12.5+2.4 

d 06:35:21.85 +03:56:27.6 4.2+1.8 15.6+2.8 

e 06:35:22.19 +03:56:29.5 2.2+1.7 10.1+2.4 

f 06:35:22.70 +03:56:33.2 2.0+1.6 4.9+2.1 

g 06:35:22.90 +03:56:28.7 15.5+3.8 52.4+5.8 

h 06:35:23.21 +03:56:17.3 1.2+1.5 3.3+1.8 

i 06:35:23.77 +03:56:24.1 1.8+1.6 5.2+1.9 

j 06:35:24.33 +03:56:21.6 4.7+2.2 20.8+3.4 

k 06:35:24.61 +03:56:20.7 2.5+1.6 10.4+2.4 

1 06:35:25.26 +03:56:21.4 2.4+1.6 8.3+2.1 

m 06:35:25.46 +03:56:21.6 1.5+1.3 4.3+1.8 

n 06:35:26.13 +03:56:20.9 2.0+1.7 8.6+2.5 

o 06:35:26.57 +03:56:23.5 1.4+1.4 3.3+1.6 

p 06:35:27.90 +03:56:18.3 3.1+1.9 11.8+2.6 

q 06:35:28.56 +03:56:25.7 3.0+2.0 10.8+2.6 

r 06:35:29.02 +03:56:14.2 12.8+3.4 42.3+5.0 

s 06:35:29.18 +03:56:18.6 2.1+1.6 6.9+2.0 

t 06:35:31.16 +03:56:15.8 1.6+1.5 4.3+1.8 

u 06:35:33.64 +03:56:16.6 3.0+1.7 11.3+2.6 

V 06:35:33.85 +03:56:18.9 2.8+1.9 7.9+2.1 

w 06:35:34.04 +03:56:20.5 2.0+1.5 5.1+2.0 

Note. — Flux is in units of 10~^* W m~^. Aperture radius 
used is 7 pixels. Flux uncertainty includes calibration uncertainty 

appears slightly deflected northward followed by a bright 
knot (g) and then a complex structure of smaller knots. 

The NIR images were flux calibrated to the 2MASS Ks 
band by determining the magnitude difference between 
the 2MASS Ks band catalog values and the ISPI image 
magnitudes for stars in common. The zero point flux in 
each NIR image after calibration is then the filter band- 
width multiplied by the 2MASS Ks-band zero point flux 
density. This results in a relation between the counts 
sec~^ in each filter image and the fiux in W cm^^. The 
narrow band continuum Kc images were scaled and sub- 
sequently subtracted from the II2 line images. Figures 2c 
and 2d show the continuum subtracted H2 1~0 S(l) 2.122 
/im and H2 1-0 S(2) 2.034 fim line images. The quality 
of the subtraction is good although there are some sub- 
traction residuals from the brightest stars present in the 
subtracted images due to differences in wavelength and 
PSF combined with changing atmospheric conditions. 

Based on the observation, this outflow has been given 
the designation MHO 1321 in the Catalogue of Molec- 
ular Hydrogen Emission-Line Objects (MHOs) in Out- 
flows from Young Star^ (D avis et al.il2009D . The fluxes 
of the individual H2 knots comprising this outflow were 
determined using a circular aperture on the continuum 
subtracted images. Table 1 lists the NIR fluxes of the H2 
emission knots. The flux uncertainty is composed of the 
rms background, poisson noise, and uncertainty from the 
flux calibration. 

3.1. IRAC color space of shocked gas 

In order to study the structure of the shocked gas, 
we applied the IRAC colo r analysis method developed 
by lYbarra &: Ladal ()2009l ). We improved the color 
analysis method by including the effects of CO 
band emission in the total emission of the shocked 

^ MHO catalogue is hosted by the Joint Astronomy Centre, 
Hawaii . http : / /www . j ach . hawai i . edu/UKIRT/MHCat / 



06:35:20.08 


+03:56:37.9 


4.1+2.0 


7.6+2.2 


06:35:21.64 


+03:56:28.8 


2.2+1.4 


8.9+2.2 


06:35:21.82 


+03:56:32.3 


3.9+1.9 


12.5+2.4 


06:35:21.85 


+03:56:27.6 


4.2+1.8 


15.6+2.8 


06:35:22.19 


+03:56:29.5 


2.2+1.7 


10.1+2.4 


06:35:22.70 


+03:56:33.2 


2.0+1.6 


4.9+2.1 


06:35:22.90 


+03:56:28.7 


15.5+3.8 


52.4+5.8 


06:35:23.21 


+03:56:17.3 


1.2+1.5 


3.3+1.8 


06:35:23.77 


+03:56:24.1 


1.8+1.6 


5.2+1.9 


06:35:24.33 


+03:56:21.6 


4.7+2.2 


20.8+3.4 


06:35:24.61 


+03:56:20.7 


2.5+1.6 


10.4+2.4 


06:35:25.26 


+03:56:21.4 


2.4+1.6 


8.3+2.1 


06:35:25.46 


+03:56:21.6 


1.5+1.3 


4.3+1.8 


06:35:26.13 


+03:56:20.9 


2.0+1.7 


8.6+2.5 


06:35:26.57 


+03:56:23.5 


1.4+1.4 


3.3+1.6 


06:35:27.90 


+03:56:18.3 


3.1+1.9 


11.8+2.6 


06:35:28.56 


+03:56:25.7 


3.0+2.0 


10.8+2.6 


06:35:29.02 


+03:56:14.2 


12.8+3.4 


42.3+5.0 


06:35:29.18 


+03:56:18.6 


2.1+1.6 


6.9+2.0 


06:35:31.16 


+03:56:15.8 


1.6+1.5 


4.3+1.8 


06:35:33.64 


+03:56:16.6 


3.0+1.7 


11.3+2.6 


06:35:33.85 


+03:56:18.9 


2.8+1.9 


7.9+2.1 


06:35:34.04 


+03:56:20.5 


2.0+1.5 


5.1+2.0 



Spitzer and near-infrared observations of a new bi-polar outflow in the RMC 



3 



01 

^ 20 
° 10 

o -10 



IRAC3.6|im " ' ^ ^ * 


V 


IRAC 4.5 ^lm 

m ■( 


^^^^^ \ 




. . 1 . . ■ . . 1 ■' . . . . 1 . . . . 1 . 


tr _ 


IRAC 8.0 urn 

1 


1 1 1 



2 10-1 210-1 

R.A. offset (arcmin) 

Figure 1. Spitzer IRAC images of the outflow. The origin is set at (a,<5)(J2000) = (06'^35™25!0, +03°56'21") 
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Figure 2. Near-infrared images of the outflow: a) H2 1-0 S(l) 2.122 /jm hnc, b) Kcont, c) continuum subtracted If2 1-0 S{1) 2.122 /xm, 
d) continuum subtracted H2 1-0 S(2) 2.034 fira. Tfie horizontal scale is in arcminutcs and the vertical scale is in arcseconds. The origin is 
set at («,<5)(J2000) = (06''35™25!0, -l-03°56'21") 



gas. Tlie distribution of the population of pure rota- 
tional levels of CO due to collisional excitation with 
H2, H, an d He was calculat e d usi ng the rate coef- 
ficients of iDraine fc Robergj (Il984f). We e mployed 
the method of iGonzalez- Alfonso et al.l ()20Q2[ ) to cal- 
culate the relative rotational population for the CO 
v=l vibrational level. The Einstein A- values for the 
CO p=l~0 rovibrational transit ions were obta i ned us - 
ing the oscillator strengths of iHure fc RoueftI (jl993D . 
In our calculations, we set uh = ?t.(H) -|- 2n(H2), 
n(He)/nH = 0.1 and n(CO)/nH = 7 x 10"^ The 
fraction of atomic to molecular hydrogen was estimated 
by considering the rate of collisional dissociation by 
H atoms, = 1-0 x 10^1° exp(-52 000/T) cm^s"! 
(|Le Bourlot et all I2002D and the rate of formation on 
grains, Rf = 3.8 x 10~^''T~^-^ cm'^s"^, derived from 



IHollenbach fc McKed (1 19 791) with the coo ling rates for 
H9 and H9O (ILe Bourlot et al.l[T999l[200l . 

Figure 3 shows the location of shocked gas in IRAC 
[3.6]-[4.5] versus [4.5]-[5.8] color space for maximum 
shock temperature of T,„ax = 6 x 10^ K for gas temper- 
atures T = 1500 - 5000 K and densities nu = 10^^ - 10^ 
cm^"^. The square brackets refer to IRAC magnitudes. 
The post-shock fraction of atomic hydrogen was found to 
be n(H)/nH 1 — 4 x lO^'^ which is co nsistent with sim- 
ulatio ns of non-dissociative C-shocks fWilgenbus et "aU 
I2000f ). In the simulations by ^ilgcnbus ct al. (2000) it 
was found that the atomic fraction is relatively constant 
over a wide range of maximum temperatures. Therefore 
we will assume our color space to be representative of 
non-dissociative C-shocks in general. The location of the 
shocked emission in IRAC color space depends on the 
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Figure 3. IRAC color-color plot indicating the region occupied 
by shocked gas composed of H2 and CO for Tmax = 6 X 10^ K 

gas density, fraction of atomic hydrogen, and the kinetic 
temperature of the gas. The [4.5]-[5.8] color is strongly 
dependent on temperature, while the [3.6]-[4.5] color has 
a strong dependence on the atomic hydrogen density. 
At high densities, the dependence on density decreases 
as the H2 gas moves toward local thermal equilibrium 
(LTE). The slope of the reddening vector is similar to 
the approximate slope of the constant temperature lines. 
Thus temperature maps of high extinction regions re- 
main accurate even if the extinction cannot be accounted 
for. However, unless extinction can be corrected for, ac- 
curate density information may not be attainable. 

We fit an analytic form to the relationship between 
color and temperature for the non-dissociative case, 

Ta = 4.19 - 0.97([3.6H4.5]) - 2.11([4.5H5-8]) 
+0.59([3.6H4-5])^ + 0.50([4.5]-[5-8])^ 

where — T/li)^ in the color space defined by 2.0 > 
;3.6]-[4.5] > -0.21([4.5]-[5.8])-h 1.5, and -0.2 < [4.5]- 
5.8] < 2.0. The difference between the analytic fit and 
the calculated temperature-color relation over the de- 
fined range is less than 10%. Additionally, one can use 
the flux in the 3.6 /im image to estimate the column 
density of the shocked H2. Using our calculations we fit 
the following analytic form to the relationship between 
column density, IRAC 3.6 /im flux density, and temper- 
ature, 

log(A^Hji^3.6) = 23.11 - 3.4OT3 + 0.742r| 



-0.0589r| 



0.071^6 + 0.012n6T3 



where N-^^ is the column density of shocked H2 in cm ^, 
F3, 6 is the IRAC 3.6 jim band flux density in units of 
MJy sr-i, and ng = nn/lO^ for 1.5 < T3 < 5.0 and 
2 < rig < 10. 

In order to investigate the color space of dissociative 
J-type shocks we simulated the gas properties of a disso- 
ciative shock by setting the maximum temperature of the 
gas to Tmax = 1 X lO'' K. Our calculations show signif- 
icant dissociation of H2 with increasing density. Figure 
4 shows the IRAC color space for dissociatively shocked 
gas. As the molecular hydrogen gets dissociated, emis- 
sion from the CO i^=l~0 band begins to dominate the 
4.5 /xm IRAC channel. 

The relationships between temperature, density and 
color are different for the case of the non-dissociative 
shock and the case of the dissociative shock. There is 
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Figure 4. IRAC color-color plot indicating the region occupied 
by dissociatively shocked gas composed of H2 and CO 

some degeneracy at the low density and low dissociation 
region of color space for the dissociative shock and at 
the high density non-dissociative shock region. These 
two models meet in a region of color space defined by 
II2 gas in LTE. Although there is degeneracy, we expect 
the distribution in color space for the dissociative shock 
to primarily lie at [4.5]-[5.8] < 0. Thus we define the 
domain of the dissociative shocked gas in color space to 
be [4.5]-[5.8] < and [3.6]-[4.5] > 1.5, whereas we define 
the color domain of non-dissociate gas to lie primarily 
at [4.5]-[5.8]> 0. By analyzing the color distribution of 
an outflow, it may be possible to distinguish between 
the cases. A pixel density distribution, produced from 
binning the colors of each pixel in the the outflow, can 
reveal the nature of the outflow by showing where most 
of the pixels lie in color space. 

It should be noted that the IRAC color analysis as- 
sumes dust and PAH emission is negligible. In order 
to prevent dust emission from affecting the color anal- 
ysis, the 8 /xm IRAC color is not used as there is evi- 
dence in some shocks of continuum dust emission within 
that wavelength r ange covered by the 8 /im channel (eg. 
iSmith et alll2006D . Additionally, PAHs are very likely 
destro yed in shocks. In a recent study, iMicelotta et al.l 
(|2010( ) show that strong shocks can destroy PAHs or 
severely denature them. 

3.2. IRAC color analysis of MHO 1321 

The IRAC 8 /im image of the outflow region reveals 
patches of absorption against the diffuse background 
(Figure 1). Of particular interest is a dark patch seen 
in absorption that bisects the outflow. We created an 
extinction map from the 8/xm data using a small scale 
median filter assuming a uniform background. We ap- 
plied this extinction map to the images of the outflow 
re gion using the m i d-in frared reddening law (KP, v5.0) 
of iChapman et all (12009). However, this is not able to 
account for the total extinction in the line of sight toward 
the outflow. Nonetheless, the temperature-color relation 
is insensitive to extinction for non-dissociative shocks. 
We estimate the background using a ring median filter 
and subsequently remove this background from the IRAC 
3.6 /im, 4.5 /im, and 5.8 /im images. A ring median filter 
is a median filter from which only the pi xels within a. n 
annulus are used in calculating the median (|Seckeiill995D . 
The scale of this filter needs to be larger than the scale 
of the shocked emission otherwise the background will be 
overestimated, yet small enough to account for the large 
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Table 2 

IRAC color analysis of H2 knots 



of the outflow knots. The distribution of pixels and the lack of 
pixels in the CO dominated region indicate non-dissociative shocks. 

scale background fluctuations. The images were shifted 
and registered with each other and then IRAC colors at 
each pixel location were determined. Figure 5 shows the 
pixel density in IRAC color space for the knots of the 
outflow. Due to the lack of pixels whose colors are in 
or near the CO dominated region and our criteria above 
for non-dissociative shocks, we conclude that this shock 
is mostly non-dissociative and we can therefore estimate 
the thermal structure based on color analysis. We com- 
pared the colors to those of non-dissociative shocked gas 
with the cutoff [4. 5] -[5. 8] < 1.5. A thermal map was 
created by estimating the gas temperature based on the 
location of the pixels in color space. Figure 6 shows the 
thermal map of the outflow. 

We find that most of the NIR H2 knots are spatially 
coincident with the high temperature regions of the flow. 
However, knots a and v do not have a corresponding 
IRAC derived temperature. The NIR images reveal stars 
in the line of sight for these knots which add to the 
emission and prevents IRAC color analysis from deriv- 
ing temperatures. Seven of the knots have estimated 
temperatures greater or equal to 3 x 10'^ K. By combin- 
ing the NIR infrared and temperature data it is pos- 
sible to estimate the extinction towards the brightest 
knots. For this we used the extinction cross-sections of 
iWeingartner Sz Draind ()2001[ ). The median extinction to 
the knots is Av = 27. The knots j and k in the vicin- 
ity of the dark clump have higher extinction compared 
to the rest of the knots. We use the median extinction 
value to de-redden the 3.6 ^m flux and use it create a 
column density map with our column density tempera- 
ture relation. Figure 7 show the column density map 
of shocked H2 in the flow. We find that there is also a 
correspondence between the NIR II2 knots and regions 
of higher column density. Using the established distance 
to the RMC of 1.6 kpc and the column density map we 
calculate the total mass of the shocked H2 (T > 2000 K) 
in the outfiow to be 3.5 x lO^o g (- 2 x IQ-^Mq). 

3.3. Outflow Source 

The source of the outflow is not seen in the NIR nor in 
the IRAC images. However, inspection of the MIPS 24 
^m image reveals a source ((a,(5)(J2000) = (06^35™25!0, 
-h03°56'21")) bisecting the outflow (Figure 8). Moreover, 
this source is spatially coincident with a dark patch seen 
in the 8 /im image. This patch is elongated nearly per- 
pendicular to the outflow and the northwest part of it has 
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3.3±0.6 


2.1±0.1 
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2.8±0.3 


2.3±0.2 
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3.3±0.8 


1.3±0.1 
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2.5±0.3 


1.1±0.2 
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2.3±0.2 


2.9±0.4 





2.3±0.1 


1.5±0.3 
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2.8±0.6 


2.0±0.2 
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2.6±0.1 


2.1±0.2 
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4.0±0.8 


1.4±0.1 
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2.8±0.6 


1.0±0.2 


t 


2.6±0.3 


1.0±0.2 


u 


3.3±0.2 


1.4±0.1 
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2.7±0.2 


2.4±0.2 



Note. — The estimated temperature is the column density av- 
eraged temperature determined from pixel colors within each aper- 
ture. The temperature uncertainty is the column density weighted 
standard deviation of the temperatures corresponding to the indi- 
vidual pixels. 

the morphology of an outflow cavity. The dark patch is 
seen in the contours that indicate mass surface densities 
obtained through extinction mapping of th e 8 /um imag- 
ing data by the method of iButler fc TanI ()2009[) . The 
contour levels correspond to mass surface densities of S 
= (2.5, 4.0, 5.0, 6.0) X lO'^ g cm'^. This smah cloud 
may be a remnant of the core from which the protostar 
formed. The morphology of the northern end of the cloud 
appears as a bi-polar outflow cavity with an opening an- 
gle e ~ 30°-40°. H2 1-0 S(l) surface brightness contours 
are overlaid on the MIPS 24 nm image that bisect the 
MIPS source and a coincident with the cavity of the dark 
cloud. This source is also detected in the MIPS 70 fim 
and MIPS 160 fim imaging data. However, we are unable 
to estimate the flux in the MIPS 70 /im band image due 
to incomplete coverage and possible contamination from 
the adjacent knot (j) whic h may have emission fro m dust 
and the 63 [O I] line (|Reipurth fc Ballvll200il ). Sim- 
ilarly, the MIPS 160 /im band image may include emis- 
sion from the knot in addition to the source. This source 
is not detected at shorter wavelengths and thus can be 
classified as a Class protostar. We propose this newly 
discovered protostar to be the source of the outflow. Ad- 
ditionally, using the column density map, we flnd that 
the mass of the shocked H2 to the east (1.8 x 10^" g) of 
this source is almost equal to the mass west of the source 
(1.7 X lO^o g) . 

Faint extended 24 fim emission is also detected at the 
locations of the brightest H2 knots (g & j). This may 
arise from fine-structur e [Fe II] lin es within the MIPS 
24 fim band (Velusam v et al.l 120071 ). This is consistent 
with our IRAC color analysis of these knots that reveal 
them to be high temperature regions (T > 3.2 x 10'^ K). 
This consistency between the IRAC color analysis and 
the MIP 24 fim emission validates our usage of the color 
space for non-dissociative shocks. 
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Figure 6. Thermal map of the outflow based on color analysis of the IRAC data. The contour levels indicate T : 
K, 4000 K. The origin is set at (a,(5)(J2000) = (06'^35™25!0, +03°56'21") 



1500 K, 2500 K, 3000 



u 
u 



o 
o 
Q 




R.A. offset (arcmin) 



-1 



Figure 7. Column density map for H2 of the outflow based on color analysis of the IRAC data. The contour levels indicate 
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Figure 8. MIPS 24 ^m image of the outflow sources. The green 
contours show the dark cloud that bisect the outflow and indi- 
cate the S values (2.5, 4.0, 5.0, and 6.0) xlO""^ g cm"^ obtained 
through extinction m apping using the IRAC 8 fim imaging data 
HButler fc Tan|[2009h . The blue contours are H2 1-0 S(l) sur- 
face brightness contours with levels indicating (0.5, 1.0, 2.0, 5.0) 
xlO~^* W m~^ sr~^. The H2 contours reveal the location of the 
outflow. The scale of the image axes is in arcseconds. The origin 
is set at (o,<5)(J2000) = (06'^35™25!0, -(-03°56'21") 



4. DISCUSSION 

4.1. Structure of the outflow 

The long axis of the flow extends to 3.3'. With a dis- 
tance of 1.6 kpe to the RMC, the flow would have a 
projected total length of 1.5 pc. The east lobe extends 
2.3' from the MIPS source to the bow shock, while the 
west lobe extends only 1'. Assuming a projected outflow 
velocity of 100 km s~^, using the east lobe we estimate 
the age of the outflow to be 10^ years. This age is consis- 



The origin is set at (a,5)(J2000) = (06'^35™25!0, -|-03°56'21") 

tent with the typical age of a Class source and thus this 
outflow may provide an accretion record of the protostar 
(IReipurth fc Ballvl[200l . We can estimate the mass flux 

of the outflow as M ~ M(ll2)vtlt where vt is the pro- 
jected outflow velcocity and k is the projected outflow 
length. Assuming the typical value of Vt = 100 km s~^ 
and using the values for the H2 mass and length of the 
east lobe, we estimate a mass flux of M ~ 10"'' Mq yi:~^ 
. This value is consistent with thos e obtained from o ther 
outflows using spectroscopic data (jPodio et al.l 120061 ). 

The NIR H2 data reveals the higher temperature re- 
gions of the outflow seen in the IRAC images. The IRAC 
images also show the cooler regions of the flow as the 
IRAC bands contain pure rotational H2 lines in addition 
to ^ = 1-0 and f = 2-1 ro-vibrational lines. 

The flow is spatially coincid ent with tw o lobes of high 
velocity CO gas observed bv iDent et al.l (|2009 ^. Simi- 
lar to many HH flows which extend further than their 
CO counterpart, we flnd the eastern half to extend be- 
yond the the eastern CO lobe. The east end of the H2 
flow ends in a large bow shock, while the west end re- 
veals a complex structure resembling either a broken up 
bow shock or multiple smaller bow shocks. Although the 
outflow is linear on large scales, the IRAC and H2 data 
reveal a region in the eastern lobe before the bow shock 
with a more chaotic structure. This deviation from a 
linear progression of knots may be due to a possible in- 
te raction with another outflow. The CO observations 
of IDent et all ([2009) reveal another flow in the NE-SW 
direction originating from the embedded cluster PL07 
(jPhelps fc Ladal[T997l ) that points toward this region. A 
collision between the two flows may explain the morphol- 
ogy and high temperature of knot r. The distribution of 
knots may also be due to variations in in jet direction 
over time and thus an indication of jet precession. 



4.2. Deflection of the Outflow 



Spitzer and near-infrared observations of a new bi-polar outflow in the RMC 
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The western end of the outflow appears slightly bent 
northward. The outflow is deflected by an angle 6^ — 20° 
where it appears to graze the densest region within the 
dark patch. The deflection angle remains small and ap- 
pears to decrease slightly beyond the interaction region. 
This is consistent with the simulations by Back et al.l 
(|2009() of outflows colliding with dense cloud cores where 
the impact parameter is large. The deflection of the out- 
flow may explain why the western end of the flow is 
shorter than the eastern end as the outflow velocity is 
expec ted to decrease after the collision (Raga & Canto 
I1995D . This is consistent with the IRAC color analysis 
that reveals high temperature shocked gas (knot j) to the 
east of the collision. 

If the outflow is composed of episodic ejection of ma- 
terial, there may be collision between clumps of mate- 
rial moving through the flow due to the velocity change 
(jRaga fc Cantd.2003 ) . As these successive clumps collide 
they may give rise to the high temperature and high col- 
umn density region (g) found slightly west of the deflec- 
tion. This interaction may also explain why the western 
lobe lacks the bow shock structure seen in the eastern 
end. 

5. CONCLUSIONS 

We present the discovery of a new bi-polar outflow 
in the Rosette Molecular Cloud and use NIR narrow- 
band and Spitzer imaging data to study the flow. We 
show that IRAC color analysis can be used to interpret 
the interaction of an outflow with its surrounding envi- 
ronment. Using our calculations of the IRAC space of 
non-dissociative shocked gas we fit analytic forms to the 
color-temperature and column density-temperature rela- 
tionships. We verify that IRAC color analysis can reveal 
regions of shocked gas and find that the NIR II2 knots 
correspond to regions of high temperature and or col- 
umn density determined through color analysis. We find 
diffuse MIPS 24 /im emission, most likely from [Fe II] 
lines, to be coincident with regions of high temperature 
thus confirming the validity of using the non-dissociative 
shock IRAC color space The NIR line ratios combined 
with the temperature estimates allow for the determina- 
tion of extinction along the line of sight which is used to 
create a column density map of the shocked II2 gas. We 
deduce that the asymmetry in the outflow is due to inter- 
actions with the dense material to the west of the outflow 
source causing deflection and possibly deceleration of the 
outflowing material. 

We thank Mike Butler for running his extinction map- 
ping code on the IRAC data and providing the mass 
surface density data shown in Figure 8. We thank the 
CTIO Blanco-4m observatory support scientists and staff 
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berto Alvarez. We thank our referee, John Bally, for his 
useful comments and suggestions. This work is based in 
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Telescope, which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided 
by an award issued by JPL/Caltech and also a NASA 
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